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Abstract 

The general decomposition of the spectral correlation tensor R^j ( k ) 
by Cambon et al. ( J . Fluid Mech. : 202 , 295; J. Fluid Mech ., 337 , 
303) into directional and polarization components is applied to the 
representation of Rij(k) by spherically averaged quantities. The de- 
composition splits the deviatoric part Hij(k) of the spherical average 
of Rij(k) into directional and polarization components H^\k) and 

II Ik). A self-consistent representation of the spectral tensor in the 
limit of weak anisotropy is constructed in terms of these spherically av- 
eraged quantities. The directional and polarization components must 
be treated independently: models that attempt the same representa- 
tion of the spectral tensor using the spherical average Hij(k) alone 
prove to be inconsistent with Navier-Stokes dynamics. In particular, 
a spectral tensor consistent with a prescribed Reynolds stress is not 
unique. 

The degree of anisotropy permitted by this theory is restricted by 
realizability requirements. Since these requirements will be less severe 
in a more accurate theory, a preliminary account is given of how to gen- 
eralize the formalism of spherical averages to higher expansion of the 
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spectral tensor. Directionality is described by a conventional expan- 
sion in spherical harmonics, but polarization requires an expansion in 
tensorial spherical harmonics generated by irreducible representations 
of the spatial rotation group SO 3 . These expansions are considered 
in more detail in the special case of axial symmetry. 


1 Introduction 

The most basic statistical property of the fluctuating velocity field in a tur- 
bulent flow is its single-time second-order correlation tensor 

Rij(x,x'-t ) = t)u!j{x\ t)). 


In homogeneous turbulence, a simpler description is possible by the second- 
order spectral tensor R tJ (k. t). which is a function of the wavevector argument 
k. Details can be found in Batchelor [1] and Craya [2], The dependence of R 1 
on the entire wavevector k and the consequent angle-dependence is involved 
in various important dynamical properties like redistribution of energy by 
the ‘rapid’ pressure-strain process. 

Many simplified models of the wavevector dependence of the correlation 
tensor R have been proposed; examples include Cambon et al. [3], Shih et al. 
[4], Ishihara et al. [5], Yoshida et al. [6], Thacker et al. [7]. These models 
share the generic form 

Rij(k,t ) = U(k,t)P i:j (k) + BU(k,t)H pq (k,t)k p k q P ij (k) 

+CU (k, t)P in (k)P jm (k)H nm (k, t ) (1) 

Complete explanation of the notation will be given later; for now, we stress 
the essential point that the anisotropic properties of the correlation are de- 
scribed by a single tensor function H(fc, t) that depends only on the wavenum- 
ber k = \k\. 

These models are revisited here by comparison with an exact decomposi- 
tion of the spectral tensor (Cambon and Jacquin [8], Cambon et al. [9, 10]) 
into terms that represent distinct properties: directional and polarization 
anisotropy. This decomposition has both a physical and a geometrical basis 
which we review. The replacement of fully anisotropic properties by spherical 

1 Here and throughout, index and index-free notation will both be used as convenient: 
R and Rij denote the same tensor, and k and k, the same vector. 
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averages necessarily restricts the description to weak anisotropy, yet even in 
this domain, Eq. (1) is inadequate because it implicitly constrains directional 
and polarization anisotropy in ways that may be inconsistent with Navier- 
Stokes dynamics. Examples are given of weakly anisotropic flows that cannot 
be described by a model of this form. Instead, we propose a new model in 
which directionality and polarization are unconstrained. The consequent de- 
scription of weak anisotropy by two tensor functions of wavenumber k is 
shown to be consistent with the dynamics. 

The paper is organized as follows. General anisotropic correlation func- 
tions are discussed in Section 2 without approximations. The fundamental 
ideas of directional and polarization anisotropy are introduced. The first 
simplification, the description of anisotropy by spherically averaged tensors, 
is introduced in Section 3. The undesirable effects of not distinguishing di- 
rectionality and polarization are described. Two examples of weak shear 
effects: sheared isotropic turbulence at short times, and the small scales in 
sheared turbulence at arbitrary times, are analyzed in Section 4, and shown 
to be consistent with a description in terms of directional and polarization 
anisotropy. Section 5 discusses the parametrization of spectral anisotropy 
by single-point moments. A refinement of the anisotropy tensor of single- 
point turbulence modeling is described, and shown to be equivalent to the 
‘structure tensor’ formalism of Kassinos and Reynolds [11], Realizability con- 
straints on anisotropic models are discussed in Section 6. They make precise 
the limitation of these models to weak anisotropy. Section 7 considers how 
the restriction to weak anisotropy can be mitigated by more accurate, higher 
order expansions based on representation theory of the rotation group SO 3 . 
Section 8 summarizes the main results. 


2 Exact relationship for arbitrary anisotropic 
second-order statistics 

The decomposition of a second-rank tensor into trace and deviator will be 
generalized to the correlation tensor, taking into account two special features: 
the solenoidal property 

kiRij(k,t ) = Rij(k,t)kj = 0 (2) 

that follows from the incompressibility of the velocity field, and the depen- 
dence on the vector argument k. In what follows, the helicity of the velocity 
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field will be assumed to vanish; accordingly, the correlation tensor is sym- 
metric: Rij(k,t ) = Rji(k,t). We recall that if the (nonhelical) correlation 
tensor is isotropic, then elementary arguments (Batchelor [1]) show that it is 
proportional to the special tensor 

R%j (kj — Sij kikj ( 3 ) 

where k = \k\ and k n = ki/k is the unit vector along k. We also recall the 
elementary property 

Pij(k)ki = P ij (k)k j = 0 (4) 

which states that P is solenoidal, and 

Pim{k)P m j{k ) = Pij{k) ( 5 ) 

which states that P is a projection. The geometric meaning of Eqs. (4) and 
( 5 ) is that at any vector fc, P(fc) is the projection onto the plane perpendicular 
to k. We will also use the obvious results 

Pmm{k) = Pmn{k)Pmn{k) = 2 ( 6 ) 

To begin, note from Eq. (2) that 0 is always an eigenvalue of R (fc,t), and 
that k itself is the corresponding eigenvector. It follows that in any frame 
centered at k in which k is one of the basis vectors, R(fc, t) can be represented 
as the matrix 

a b 0 

R= b c 0 (7) 

_ 0 0 0 _ 

characterized by exactly three real scalars, where we again remark that the 
absence of helicity implies that R is symmetric in any basis. A trace-deviator 
decomposition in the plane normal to k yields 

e 0 0 d 6 0 

R= 0 e 0 + 6 -d 0 (8) 

_ 0 0 0 J [ 0 0 0 _ 

where 

e = |(a + c) d = |(a — c) (9) 
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In this frame, P, as the projection onto the plane perpendicular to k is 
represented by the matrix 

'10 0 ' 

P= 0 10 (10) 

_ 0 0 0 _ 

Introducing the two independent symmetric matrices 

'1 0 0 1 [010' 

Mi = 0-10 M 2 = 10 0 (11) 

_ 0 0 0 J [ 0 0 0 _ 

we therefore have 

R = eP + dMi + 6M 2 (12) 

In terms of the complex quantities 

Z = d + ib M = Mi — iM 2 (13) 

Eq. (12) can be written alternatively as 

R = eP + 5ft (ZM) (14) 

An obvious coordinate system in which k is a basis vector at every k is 
the spherical coordinate system in fc-space, in this context called the Craya- 
Herring frame (Craya [2], Herring [12]). For the application of this coordinate 
system to explicit expressions for Mi and M 2 in terms of the helical mode de- 
composition using the (complex) eigenvectors of rotations about k, see Cam- 
bon and Jacquin [8] and Waleffe [13]. Eqs. (12) and (14) express R in terms of 
the minimal number of scalars: the three real quantities b(k. t). d(k , t), e(fc, t) 
or equivalently, e(k,t ) and the complex scalar Z(k,t). 

Defining 

d, 6 0 

R po * = dMi + 6M 2 = 6 -d 0 (15) 

_ 0 0 0 _ 

the decomposition in Eq. (12), 

R = eP + R po * (16) 
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is characterized by the properties 

e = IR : P = |tr R R pul : P = 0 (17) 

Because these properties are independent of the coordinate system, we can 
also arrive at Eq. (16) by coordinate-free arguments. Thus, define the pro- 
jection of R along P by 

R P = i(R : P)P (18) 

where the factor of 1/2 is due to Eq. (6). The operation so defined is a 
projection because [R p ] p = R p . Accordingly, the decomposition 

R=i(R:P)p+[R-i(R;P)p] (19) 

coincides with Eq. (16) after introducing the definition Eq. (17) of e and 
replacing Eq. (15) by the coordinate-free definition 

R po1 = R — |(R ; P)P (20) 

The polarization tensor is geometrically very simple: R po1 has one zero 
eigenvalue because it is solenoidal, and since R and R po1 are both solenoidal, 

tr R po1 = R po1 : P = 0 (21) 

Accordingly, its characteristic polynomial is simply 

p( A) = A 3 - (iR po/ : R po/ )A 

which also follows directly from the explicit expression Eq. (15). It follows 
that the eigenvalues of R po1 are 0, R po1 : R p ° l . Since the eigenvalues of eP 

are obviously just e, e, 0, the eigenvalues of R are e ± \J\R po1 : R po/ , 0. The 
realizability of R is therefore simply the condition 

e > \/\R po1 : R po1 (22) 

Note from Eqs. (13) and (15) that 

ZZ* = \R po1 : R pul (23) 

Thus, although the scalars d = 5RZ and b = ^sZ are coordinate-dependent, 
the magnitude |Zj is a geometric invariant. 


6 



Eq. (19) is a straightforward generalization of the trace-deviator de- 
composition in which the trace, the projection along 5 t j. is replaced by the 
operation of projection along P. 

The general decomposition Eq. (16) can be rewritten in a form that 
isolates the purely isotropic part by projecting e onto its spherical average, 

f/ ( /c,t ) = 4 ^£ fc e(M)d2fe (24) 


where $ {...)d 2 k denotes integration over a spherical shell of radius k. Note 

J Sb 

trivially U(k,t ) = - — — <1 U(k,t ) d 2 fc, Eq. (24) does define a 

47T K J Su 


ls k 

that since 


projection. Defining 
and 

we have 


£(k,t ) = e(fe,t) — U(k,t ) 
R dir (k, t) = £(k,t)P(k) 


(25) 

(26) 


R(M) = U(k,t)P(k) + R dir (k,t) + R t) 

Isotropic part Directional anisotropy Polarization anisotropy 

(27) 

The decomposition of the correlation tensor in Eq. (27) has a simple but 
important geometrical significance: just as in the trace-deviator decomposi- 
tion, each component transforms into another term of the like form under 
arbitrary rotation of axes. 


3 Description by spherical averages alone 

Even in single-point turbulence modeling, it is often useful to introduce 
an ansatz for the energy spectrum. For example, such assumed spectra 
have proven useful in closing the ‘rapid’ pressure-strain correlation [3, 4], 
But the obvious practical difficulty of a complete description of spectral 
anisotropy has motivated a search for simpler descriptions of anisotropy. 
The approach adopted here is to reconstruct the angular dependence of the 
spectrum through appropriate tensorial expansions, using either functions of 
the wavenumber alone, as in the expression R (k,t) = U(k,t)F(k,H(k,t)), 
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or even more simply, functions of single-point moments as in the expression 
R (k,t) = U(k, t)J-(k, b (t)): in each case, E denotes a dimensionless isotropic 
tensorial function of its arguments. The presence of k in the list of argu- 
ments distinguishes these expansions from the tensorial expansions familiar 
in single-point modeling. The Introduction briefly described some models of 
this type. This type of modeling will be reconsidered in this section in the 
light of the distinction between directional and polarization anisotropy. 

By spherical integration, we can construct two obvious tensor measures 
of anisotropy that depend only on k: they are defined by 


2E(k,t)H^(k,t) 

2£(M)H (z) (M) 

where 

E(k,t) = e(k, t)d 2 k = 4TYk 2 U(k) (30) 

is the energy spectrum. 2 The notation follows Cambon and Jacquin [8], and 
is motivated by the characterization in Section 2 of directional anisotropy 
by the scalar e(k,t ) and of polarization anisotropy by the complex scalar 
Z(k,t). Obviously, 


R dir (k,t)d 2 k 

■kSk 

f R pol {k,t)d 2 k 

J Sb 


(28) 

(29) 


2E(k,t)tv H ( e \k,t) = 


tr R dtr (k,t)d 2 k = (f (e-U) tr P (k)d 2 k = 0 (31) 
ks k 


and, in view of Eq. (21), 


2E(k,t)ti H {z) (k,t) = l tr R pol (k,t)d 2 k = 0 (32) 

Js k 


so that both of H^ e,z ^ are trace-free. 

We wish to construct a modeled correlation tensor that depends only 
on the spherical averages Y\^ e,z \ The discussion in Section 2 motivates con- 
structing R dtr and R po1 separately. To begin, note that R dlT depends linearly 
on H( e ) and is proportional to P. The simplest assumption consistent with 
with these properties is 

R dir (k,t) = AU(k,t ) [H (e) (M) : P(Jfe)] P(fe) (33) 

2 Our definition of U is somewhat nonstandard; it is customary to define it so that 
E = 2 t xk 2 U. 
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with an undetermined constant A. Equivalently, in terms of the solenoidal 
tensor PH^P, Eq. (33) sets R dir = A\yY\^P] dir . The constant A should be 
chosen to be consistent with the definition Eq. (28); the spherical average 3 
of each side of Eq. (33) gives 

2£(M)H (e) (M) = ~^AE(k,t)H^(k,t) (34) 

so that A = —15. 

The treatment of the polarization tensor is somewhat less straightforward. 
R po1 must be solenoidal and linear in H^. These requirements suggest the 
form R po/ = PH^P. But in addition, we must take into account that R po/ : 
P = 0. A general form consistent with all constraints is therefore 

fig'(M) = BU(k,t) t)p„j(k) - iH $(k,t)P„(k)P tj (k) 

(35) 

with an undetermined constant B. Again, in terms of the solenoidal tensor 
Phb 2 )P, Eq. (35) sets R po1 = B[PY\^P] po1 . Spherical averaging as in Eq. (34) 
gives 

2E(k,t)HW(k,t) = ±BE(k,t)HW(k,t) (36) 

so that B = 5. 

Combining the results of Eqs. (33) and (35), we obtain the required 
representation 

R ij (kJ) = U(kJ)P ij (k) + 15U(k,t)P ij (k)H^{k,t)P„(k) 

+5 U(k,t) [P irx (k)P jm (k)H^(k, t) - iP mn (k)HM(k,t)P tj (k)} (37) 

or equivalently, 

R.jtkJ) = U(k,t)Pij{k) - 1 5U(k,t)P ij (k)H^(k,t)k r k q 

+5f J(k,t) [P tn (k)P jm (k)H^(k,t) + i P«(fc)tfW(M)M,] (38) 

This equation is the main result of this paper. It shows that a completely con- 
sistent description of weak anisotropy without arbitrary constants is possible 
using independent descriptors of directionality and polarization by tensor 
functions of the wavenumber k alone. 


3 The derivation requires the formulas (Cambon et al. [3]) 




kikjd 2 k = and 


kikjk m k n d k — , \dpj d lnn dimdnj dj n d rn j\ . 
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Note that since tr bfi e ) = 0, the quantity H^(k,t)k p k q that appears in 
Eq. (38) is a second order spherical harmonic: after choosing a polar axis and 
introducing spherical coordinates, it would be expressed in terms of Legendre 
functions in the standard way. The term in brackets containing has an 
analogous interpretation as a Censor spherical harmonic (Zemach, [14]). We 
will follow Zemach [14] in referring to scalar spherical harmonics (SSH) and 
tensor spherical harmonics (TSH) henceforth. From this viewpoint, Eq. (38) 
states the lowest order terms in expansions of R dzr and R po/ respectively in 
scalar and tensor harmonics, or in irreducible representations of the rotation 
group SO 3 . This connection will be developed further in Section 7. 

Eqs. (28) and (29) give expressions for the spherically averaged correla- 
tion function 

<, Pij{k,t ) = jf Rij(k,t)d 2 k = 2 E(k,t) + H\f(k,t) + H[f(k,t)} . (39) 

Eq. (38) states the important conclusion that a solenoidal tensor R can 
indeed be constructed from the spherical averages this conclusion is 

not obvious because the solenoidal property is lost on spherical averaging. In 
modeling, it has been the general practice to characterize anisotropy by one 
tensor instead of two. Thus, many models are based on 

(Pij(k, t) = Rij(k , t)d 2 k = 2 E(k, t ) U5 tJ + H i:j (k, t )) (40) 


where H = + H^. Such modeling mixes directional and polarization 

anisotropy. We next ask whether a solenoidal correlation function can be 
constructed consistent with H alone. 

Proceding as before, we set 


p .dir 




AU{k,t)(H : P)P 
BU(k,t ) fPHP 


i(H : P)P] 


(41) 


We already note that these equations implicitly impose some relation between 
R d * r and R poi and therefore cannot be entirely satisfactory. On spherical 
averaging, we find 

H = (A- 4 + ^) H < 42 ) 

The solution is not unique: it is 


A = 15 + 3a B = —a 


(43) 
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where a can be a function a = a(k,t). Thus, the spectral tensor is 

R (k,t) = U(k,t)P(k) + (1 + ±a)U(k,t)(H(k,t) : P(fe)P(fe) 

-\aU{ k, t ) [P(jfe)H(fc, t)P(k) - ±(H(Jfe, t) : P(fe))P(fe)] (44) 

and 

H (e) = (l + | a) H, H^ = -faH (45) 

which implies the relation 

faHO = -(1 + §a)HW. (46) 

Eq. (44) is exactly the proposal of Cambon et al. [3] for the first general 
model linking R with ip. The aim was to derive an equation for the spherical 
average <p(k,t) from the evolution of R in the presence of mean flows with 
constant mean velocity gradients. The basic problem is that spherical aver- 
aging introduces a closure problem for the pressure-strain and transfer terms, 
so that a model for R(fc) was needed, parametrized by functions of k only, 
and so related in a straightforward way to ip(k, t). Starting from an equation 
for R (k,t) that includes the exact linear terms (recalled at the beginning 
of the next section) and contributions from triple correlations closed by an 
anisotropic EDQNM theory, Eq. (44) allowed Cambon et al. [3] to system- 
atically derive a closed equation for (p(k, t ), but involving a(k, t) from (44) as 
an ajustable parameter. Other more empirical models in which anisotropy is 
parametrized entirely by a tensor function of wavenumber k have been pro- 
posed [15, 16]; they explicitly use some approaches of single-point turbulence 
models to close the pressure-strain correlation. 


4 Application: short and long time behavior 
of weakly sheared turbulence 

We have noted that Eq. (38) can be understood as a lowest order expansion in 
anisotropy valid only for small departures from isotropy; the restrictions will 
be clarified shortly by analysis of the realizability of the correlation proposed 
in Eq. (38). This section will give examples of weakly anisotropic flows that 
permit direct analysis, and will show that the correlation tensor is indeed 
described by Eq. (38). The first example is homogeneous shear flow with 


11 



isotropic initial conditions treated by rapid distortion theory (RDT) in the 
short-time limit in which anisotropic effects are necessarily small. 

Using the notation of Cambon and Scott [17], we begin with the general 
RDT equation 

f)R dk r) 

= -^-—R t] (k,t) - M m (k)R nj (k,t ) - M jn (k)R in (kR) (47) 

where = (Si m — 2kik m )A m j, A im = dUi/dx m is the mean velocity gradient, 
and the wavevector k n satisfies dkij dt = —Ajikj. We will consider evolution 
away from an isotropic initial condition R(fc, 0) = t/(fc)P(fc); we analyze the 
evolution at very short times, when the effects of shear remain weak. 

Consider the first order Taylor series expansion R(fc, t ) = U(k)P+tR(k , 0). 
Evaluating Eq. (47) at t = 0 using the isotropic initial condition R(fc,0) 
= U(k)P(k) leads easily to 

■df {M) = A pnKknkU\k)P ZJ (k) + U{k)k i k p A pn P nj (k) 

+U (k)kjk p A pn P ni (k) - U(k)A in P nj (k ) - U(k)A jn P ni (k ) 

= S pn k p k n kU'(k)P tJ (k) - U{k)P ip {k)P nj {k)S pn (48) 

where S tJ = \{Aij + Aji) is the strain rate. It follows that 

Rf(k, t ) = \tS pq k p k q [kU'(k) + U(k )] P tj {k) (49) 

H l lt(k, t ) = -tU(k) [P im {k)P jn (k)S mn + \S pq k p k q P l3 {k)\ (50) 

We remark that the U' term comes from the d/dk effect in the RDT equa- 
tions: it is a conservative linear energy transfer mechanism in fc-space, which 
therefore appears as a directionality effect; polarization effects instead arise 
when energy is transferred between different tensor components of the cor- 
relation. We see from Eqs. (49)-(50) that both effects are relevant. 
Spherical integration gives 

= + and = (51) 

This leads to a(k ) = 15/(5 + (k/ E)dE / dk) in Eq. (45), which is constant in 
any fc-space region where E obeys power-law scaling. 
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It is useful to supplement this short-time analysis by long-time nonlinear 
analysis, as proposed by [18, 19, 20] and others. In these computations, it is 
assumed that at least for some range of scales, the shear may be treated as a 
weak perturbation of isotropic turbulence: the calculation itself will identify 
an appropriate small parameter. The most recent such analysis, by Yoshida 
et al. [5, 6] overcomes some limitations of previous work, and introduces a 
Lagrangian viewpoint, with important conceptual and computational advan- 
tages. Yoshida et al. refer to this calculation as ‘linear response theory’ for 
turbulence, since it is accomplished by linearizing about an isotropic nonlin- 
ear state. As in the RDT problem above, anisotropy is therefore weak, and 
we can again expect Eq. (38) to describe the spectral tensor. 

General kinematic considerations (the same as in [3, 4]) lead again to 

Rij(k) = U{k)P ij {k) + b(k)P ij (k)S pq k p k q + 2a\k)P in (k)P jm (k)S nm . (52) 

from which we immediately deduce 

= (b — a')(S pq kpk q )Pij(k) 

R P ij l = 2 a' ^P irn (k)P jn (k)S mn + lS pq k p k q P ij (k)'j (53) 

exactly the same structure as in Eqs. (49)-(50) except for the scalar functions 
of k ; consequently, 

Hif = H\f = -I WSij (54) 

where O^ik) = [b(k) — a'(k)]/U(k) and 6^ z \k) = 2 a'(k)/U(k) are time scales. 
In a Kolmogorov inertial range, we will have 

0 (e) = (A - B)e~ 1/3 k ~ 2/ 3 9 {z) = Ae- 1/3 k~ 2/3 (55) 

where A and B are universal constants. In this case, 

|AH^ = A(5_ a)H^ (56) 

and the anisotropic part of the spectrum satisfies 

£H (e) ~ £H (z) ~ k~ 7/3 (57) 

the scaling suggested by dimensional analysis since H( e ) oc S (see Yoshida et 
al. [6] for more details and references to earlier work). 
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Yoshida et al. find theoretical values from spectral closure theory A « 
— 0.16 and B & —0.40. Experimental and DNS measurements give values 
closer to A ps —0.12 and B ~ —0.009. It is found that making realistic 
corrections to the theoretical values to account for the finite inertial range 
in the measurements, results in much closer agreement. However, for our 
purposes, the actual values are not so crucial; the important observation is 
that a weakly anisotropic spectral model based on a single tensor H imposes 
some fixed proportionality = AH^ as shown by Eq. (46); such models 
cannot be consistent with both the short- and the long-time results Eqs. (51) 
and (56). On the other hand, the model Eq. (38) can be consistent with 
both limits. 

In the case of pure rotation, where A is antisymmetric and hence the strain 
S vanishes, no information is given by the previous ‘linear’ approaches. Nev- 
ertheless, rotating turbulence provides another example of how directionality 
and polarization must be separated in general. The inviscid RDT equation 
(e.g. Cambon and Jacquin [8]) for this problem implies 

R dir (k,t) = R dir (M), R pol {k,t) = exp(4iflk 3 t)R pol (k,0) (58) 
which yields 

= Hff(kA)-* 0. (59) 

The damping of polarization anisotropy reflects the angular phase mixing 
due to the anisotropic dispersivity of inertial waves. We see then that the 
kinematics of turbulence under rapid rotation is dominated by directional 
anisotropy alone. Spherically averaged polarization can be neglected even 
in the presence of nonlinearity, but directional anisotropy can be created by 
nonlinearity, even if it is initially zero, in an incomplete transition from 3D 
to 2D structure [8, 10, 21]. 

5 Parametrization by single-point moments 

In Reynolds stress modeling, the basic descriptor is the anisotropy tensor b 
which can be recovered in the present notation by 

poo r oo 

/ dk b(t)E(k,t) = / dk E(k,t)H(k,t) (60) 

Jo Jo 
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with the classical definitions for kinetic energy /C and Reynolds stress anisotropy 
tensor b: 


I POO j 

*#) = = / E (k, t)dk, bijit) = ==f - -5ij 

Z •'O a n a n 0 


An example of a parametrization of the correlation tensor by single-point 
moments is the proposal of Shih et al. [4]; written in terms of directional 
and polarization anisotropy components, it is 



Rijik) = U(k,t)P i:j (k) + ip 3 (k,t)P ij (k)('y - \)b pq k p k q 
ik,t)Pi n ik)Pj m ik) (b nm T 2^nmb pq k p k q ') . 

(61) 

This expression can be identified with Eq. (38) with the particular choices 


u 2> E(k, t) E(k,t ) 

(62) 

so that 

5H (e) = 3(7- i)HW 

(63) 


Since this result fixes a definite proportionality between H ^ and H^, how- 
ever, it is really a special case of Eq. (44), so that the constant 7 could be 
related to the parameter a. It is not necessary to give the explicit relation; 
the important fact is such models cannot be consistent with the equations 
of motion or the general kinematics. Considerations of realizability which 
will be discussed later, led to the choice 7 = |, which removes directional 
anisotropy entirely, which is clearly inconsistent with the problems of shear 
and rotating turbulence just analyzed. The important conclusion is that the 
spectral tensor cannot be uniquely reconstructed from the stress anisotropy 
alone. 

Another model of this type recently proposed by Thacker et al. [7] has 
the form 


Rij(k,t ) = U(k,t)Pij(k ) + P in {k)P jm (k)U a (k,t)[b nm (t) + \b nm b pq (t)k p k q \. 

(64) 

This model is essentially the same as Eq. (61) with the special choice 7 = |, 
which removes directional anisotropy. 
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The fundamental decomposition of anisotropy into directionality and po- 
larization suggests introducing two additional single-point descriptors 



dk b (e) (t)E(k,t) 
dk b {z \t)E(k,t) 



dk £(M)H (e) (M) 
dk £(M)H w (M) 


(65) 


so that 

b = b (e) + b w (66) 

which can be considered a refined decomposition of the Reynolds stress 
anisotropy. 

As an example of Eq. (66), consider the short time RDT results of the 
previous section. Multiplying each result in Eq. (51) by 2 E and integrating 
over k gives 

z,( e J _ _2_ c..f u( z ) _ _ 2 c , 

~ 15 ^A ij ~ 15 ^A ij ~ 5WA 

so that b^ = — (1/3) b^^, or equivalently, — 2b^ = b. This condition, re- 
ferred to by Kassinos et al. [11] as ‘dimensionality and componentality having 
the same anisotropy’ although strictly derived only at short times, is con- 
sidered to be valid at large time provided that the mean flow is irrotational. 
Comparison with Eq. (51) shows that Eq. (67) certainly does not hold for 
spectral quantities, because the ratio of hfl e ) to H^) at the first order in time 
depends on the initial spectrum through the ‘linear transfer’ generated by 
the time dependence k(t). 

We can also compute b*- 6 ’^ for the results of Yoshida et al. [6]. To obtain 
a definite result, we will follow previous calculations of this type ([18, 19, 
20] to name just a few) and integrate the inertial range spectrum E(k ) = 
C^e 2 /^ -5 / 3 over a range of scales k > fco, where A;^ 1 is an integral scale of 
the turbulence. The result is 



iff = ~ A)e-' , % 2/3 S, j 

iff = (68) 

In comparing the results Eqs. (67) and (68), we see that the time-scale in 
short-time RDT is simply elapsed time t, whereas it is a turbulent time- 
scale oc e -1 / 3 A;o 2 ^ 3 in the linear response theory. Moreover, the ratios of 
components b\f /bif are not the same. 

We would like to add some brief remarks on the use of the refined single- 
point anisotropy measures in turbulence modeling. One attempt to 
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improve the prediction of RANS models, and especially to enforce greater 
consistency with rapid distortion theory, has been to introduce new ten- 
sors in the models. After Cambon et al. [9], the most prominent exam- 
ple is the structure tensor formalism proposed in Kassinos et al. [11]. Two 
kinds of anisotropy, called ‘componentality’ and ‘dimensionality’ were consid- 
ered. Since ‘componentality’ is nothing else than the Reynolds stress tensor 
anisotropy identified by b tJ . new information is carried only by the ‘dimen- 
sionality’ tensor 

Dij = 2 J e(k 1 t)k i kjd 3 k. 

A trace-deviator decomposition for this tensor can be expressed as 

Dij = 2/C - 2 4 e) ) , (69) 

so that the ‘dimensionality’, measured by the deviatoric part of D^, is in- 
duced by directional anisotropy alone, and ‘componentality’ mixes direction- 
ality and polarization, via b = + b^. The reader is referred to Salhi and 

Cambon [22] for more complete relationship with ‘circulicity’ and ‘strophol- 
ysis’ tensors from Kassinos et al. [11]. 


6 Realizability 


The construction of a modeled correlation tensor in terms of spherical aver- 
ages leaves open the question whether the result Eq. (38) in fact describes a 
possible correlation. This is the issue of realizability. Obviously, in the con- 
text of weak anisotropy, the anisotropic part is a small perturbation of the 
isotropic part [7(/c, f)P(fe), and as this part is realizable, it remains realizable 
under sufficiently small perturbations. 

Realizability is imposed in RANS models through the positivity of the 
Reynolds stress. But in a model in which directionality and polarization are 
treated separately, it is natural to impose realizability conditions on each 
component. Recall in the case of Eq. (38), that e — U (l — 15 H^k p k^j. 
Thus the simple necessary realizability condition e > 0 is violated if the 
largest positive eigenvalue of H.\f is larger than 1/15: this tensor being 
tracefree, at least one positive eigenvalue must exist; 4 therefore, this realiz- 


Tn fact, from 


e(k , t)Pij(k)d 2 k 


= 215(A) ( - 6 ., 


+ H ^ 

n ' n ij 


these eigenvalues are 


17 



ability constraint quantifies precisely how ‘small’ the anisotropy must be to 
admit description by Eq. (38). 

These difficulties are not unexpected, because we have formulated an 
expression applicable only to weak anisotropy. The most direct remedy is 
to introduce higher order angular harmonics in the expression for R. This 
problem is considered in the next section. We thus find that strong anisotropy 
is reflected not only in large values of the lowest order harmonics but also by a 
large number of angular harmonics expansion needed for a correct description 
of e — Z angle- distribution. 

So far, we have only treated the condition e > 0; the stronger condition 
e> \Z\ is not discussed here for the sake of brevity. We can just note that the 
condition e > 0 which depends only on the magnitude of H.\f in (38) is the 
most sensitive condition; this could explain why the directional anisotropy 
is neglected in some ad hoc models, for examply by choosing 7 = 1/2 in Eq. 
(61): even if is never zero in ‘true’ anisotropic homogeneous flows, it 
is safer to suppress its contribution in some spectrum models to avoid too 
much sensitivity to realizability issues. Even the value of the constant in the 
quasi-isotropic model by Launder et al. (1975) [23] seems (fortuitously ?) to 
reflect this condition, as shown by [3, 4], 

This problem does not occur for simpler models based only on H or b, both 
of which can be computed for any anisotropy whatsoever, although there is 
no guarantee that the predicted dynamics will be realistic (Rubinstein and 
Girimaji [24]). Thus, the inclusion of will improve a model only for weak 
anisotropy, as the low order suggests. Although, as noted earlier, it is possible 
to derive from Eq. (38) closed equations without adjustable parameter for 
E(k,t ), H( e )(/c,t), H ( z \k,t) from any closed equation for R tJ . and the same 
procedure also yields closed equations for single-point moments /C, b^, b^, 
the disappointing result, found independently by Cambon et al. [9] and 
Kassinos et al. [11], is that the resulting model behaves in the RDT limit 
worse than a conventional RST model, because of possible loss of realizability 
when the anisotropy is not necessarily small. 

bounded by ±1/3, as they are for any deviatoric tensor derived from a definite-positive 
matrix. 
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7 SO 3 decomposition 

The realizability constraint restricts the application of a theory based on 
spherical averages alone; this section will describe how to construct more 
accurate approximations using higher order expansions based on irreducible 
representations of the rotation group SO 3 . This representation theory has 
recently been found very useful in clarifying the scaling properties of corre- 
lation functions in turbulent flows (Arad et al, [25]). 

7.1 Directional anisotropy 

The spectrum models that have been considered so far suggest how the 
anisotropic part of the correlation can be expanded in powers of k. For 
example, we can continue the development of the first term on the right side 
of Eq. (38) beyond the second order by writing 

£ = U^ n (k, t)k m k n + U^ nrs (k, t)k m k n k r k s + ... (70) 

where U 2 = — 15H( e ). The expansion is restricted to polynomials of even 
order because of the parity property Rij(-k) = Rij(k). 

Although the notation suggests that Eq. (70) proceeds in powers of hi, 
it does not do so without some restrictions on the coefficients. For exam- 
ple, U^ nnrs — 5 mn A rs + • • • will generate a term A rs k r k s that could be in- 
cluded in the U 2 term. Without presenting any details, suffice it to say 
that if all such redundancies are eliminated, then U 4 will belong to a 9- 
dimensional representation of SO 3 on homogeneous quartic polynomials sat- 
isfying y 2 Ufj mn kikjk m k n = 0. This is discussed in detail in standard ref- 
erences such as Weyl [26]. Equivalently, after choosing a polar axis n, the 
expansion could be described in terms of Legendre functions as a spherical 
harmonics decomposition following Cambon and Teissedre [27], 

No n 

£( fe ) = E E eZ{k,t) P^icos 9) exp{xm(p)' 

n=lm=—n v 

YSilfiw) 

where 9 = arccos(n*fc) is the polar angle and <p is the azimuthal angle in a 
system of polar-spherical coordinates with axis n, and P 2 ™ are the associated 
Legendre polynomials of degree 2 n and order m. 


(71) 
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7.2 Polarization tensor 


The analogous higher order expansion of the second term on the right side 
of Eq. (38) is 

R% l (k,t) = § [P lm (k)P jn (k) + Pin{k)P jm (k) - Pij(k)P mn (k )] x 

P Tmnrs(k, t)k r k s + ' ' '} (72) 

where T° = 5hb 2 b 

The structure of the higher order terms in the expansion of polarization 
is not so simple as the higher order terms in the expansion of directional 
anisotropy, and the steps which make Eq. (72) an orthogonal expansion in 
irreducible representations are less obvious and standard than the steps lead- 
ing from Eq. (70) to Eq. (71). Although T° (oc hb 2 )) is simply a constant 
trace- free second-rank tensor, T 2 consists of matrices with quadratic polyno- 
mial entries; their decomposition into irreducible representations or tensor 
spherical harmonics (TSH) can be summarized as follows (compare also the 
discussion in Arad et al. [25] for more general tensor quantities): invariant 
tensors of differential operators can be constructed. Their action on homoge- 
neous polynomials belonging to a representation on scalar functions generate 
the appropriate representations on tensors. In the case just mentioned, we 
find that T 2 belongs to a 25-dimensional representation that is decomposed 
into irreductible representations of dimensions 1,3, 5, 7, 9. The representation 
of dimension 1 cannot contribute to polarization, and the representation of 
dimension 3 does not survive solenoidal projection. The remaining repre- 
sentations of dimensions 5,7,9 define solenoidal TSH. In the interest of con- 
creteness, we list the TSH that correspond to the irreducible representation 
of dimension 5. To compute it, we use the infinitesimal generators of SO 3 : 

L x = kydjdk z — k z d/dk y 

L y = k z d/dk x — k x d/dk z 

L z = k x d/dk y — k y d/dk x (73) 

Let P -2 denote any five harmonic quadratics, for example, {k x — k y , A; 2 — 
A; 2 , k x k yi k y k z , k z k x } The representation of dimension 5 is defined by 

Nh = ( 'LiLj + LjLi — | L p L p 5ij ) P/ 7 = 1, ...5 (74) 
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Examples are the matrix so generated from — ky\ 


m - K) 

^k x ky 

3k x k z 


2k x k y 
m - ki ) 

3k y k z 


3k x k z 
3ky k z 




K)\ 


( 75 ) 


and two others obtained by simultaneous cyclic permutations of x,y,z and 
of tensor indices (note that the three tensors so generated sum to zero), and 
the matrix generated from k x k y 


2 k x k y 

K + ky- 2kl 
3k y k z 


h 2 - 4 - h 2 — Oh' 2 
^ x ' ™ y ^ z 

2k x ky 

3k x k z 


3 kyk z 

3 k x k z 

4 k x ky 


(76) 


and two others obtained by the same cyclic permutations. 

This analysis shows that it is only a coincidence that at the lowest order of 
anisotropy, directionality and polarization are both described by a constant 
trace-free second-rank tensor. At even the next order, the descriptions are 
quite different: as noted above, directionality is described by a homogeneous 
fourth degree polynomial, but polarization is described by three distinct types 
of second rank tensors with quadratic polynomial entries. 

We believe that these expansions in scalar and tensor spherical harmonics, 
which clearly have very different mathematical origins, again underscore the 
difficulties of confounding directional and polarization anisotropy along the 
lines of the standard anisotropic models. 


7.3 Axisymmetric turbulence 

The problems raised by extending the description from weak anisotropy to 
arbitrary anisotropy can be better understood in the comparatively simple 
case of axial symmetry. Let the unit vector n be the axis of symmetry. In this 
case, any tracefree tensor function of k alone obtained by spherical averaging, 
or any single-point moment obtained by integration over all Fourier modes 
can be expressed as Hij = | H n (3nirij — 5ij ) in terms of the single axial 
component H n = HijUiUj. As for the fc-dependent spectra, a polar-spherical 
system of coordinates can be introduced, so that Eq. (71) reduces to 

Vo 

£(M) = £ e 2 n (k,t)P% n (cos9) (77) 

71 = 1 
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in which only terms with m = 0 appear. In addition, the polarization 
anisotropy also can be expanded as 

Ny 

Z(k,t) = sin 2 9 ^ z n (k, i)P°(cos 9) (78) 

n = 0 

in which the complex- valued Z is defined as in Cambon and Jacquin [8], us- 
ing the symmetry axis n to define the helical modes. In Eq. (78), z n is real 
for even n and purely imaginary for odd n due to the hermitian symmetry 
property Z(—k ) = Z*(k ); imaginary terms can represent the breaking of 
mirror symetry present, for example, in rotating turbulence and also yield 
some ‘stropholysis’ terms (Kassinos et al. [11]). If axial symmetry is un- 
derstood to include invariance under reflections through planes containing 
the polar vector n, then Z is real, and only terms of even order appear in 
the expansion of Z. Restricting to even orders with Ny = 21V 0 , Eqs. (77) 
and (78) can be recovered from Sreenivasan and Narasimha [28] and from 
Cambon & Teissedre [27], 

Note that the factor sin 2 9 is essential in the expansion of Z in Eq. (78) 
because polarization anisotropy must vanish when the wave-vector is parallel 
to the axial vector n : reference to Eq. (15) shows that R po1 can only be 
axisymmetric if b = d = 0. The expansion of Z in Eq. (78) therefore is 
somewhat special. The general spherical harmonics expansion for Z in [12] 
is not consistent with this property. At the lowest order (7V 0 = 1, Ny = 0), 
H-^kikj = with P 2 ° = 1(3 cos 2 9 — 1), and Z = | sin 2 9 1 defining 

Z as in Cambon and Jacquin [8]. 

Since the spherical harmonic decomposition required for the expansion of 
directionality to higher order is entirely standard, we consider the expansion 
of polarization to higher order. We will construct the axisymmetric tensor 
spherical harmonics of the first order beyond HC), without however present- 
ing derivations. We noted at the end of Section 7.2 the existence of TSH 
belonging to the irreducible representations of dimensions 5,7, and 9. We 
find exactly one axisymmetric TSH belonging to each of these representa- 
tions. That corresponding to the irreducible representation of dimension 5 is 
generated by the (unique) axisymmetric quadratic polynomial 2 k 2 z — A; 2 — k 2 
by the differentiation process described in Section 7.2. The corresponding 
TSH is 

—2k'l + 4 ky — 2 k 2 z —6k x k y 3 k x k z 

A 5 = —6k x k y 4k 2 x — 2k 2 — 2 k 2 3 k y k z (79) 

3 k x k z 3 k y k z —2 k 2 — 2 k 2 + 4k 2 
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so that 


( 80 ) 


A 5 : P = 2 k~\(kl + *J) - 2 k\) 

Recalling the relation Eq. (23), 

Z = A :~ 4 ^6(k 2 + ky ) 2 + 6(A c 2 + ky)kfj = 6 sin 2 9 (81) 

in agreement with (78) for n = 0. 

The TSH belonging to the irreducible representation of dimension 7 prove 
to change sign under inversion through planes containing the polar axis. We 
simply note the (again unique) axisymmetric TSH 

2 k x k y -k'l + A; 2 -Ak y k z 

A 7 = —A; 2 + ky —2k x k y 4 k x k z (82) 

—Ak y k z ^k x k z 0 

for which e = 0 and Z = lChA; -3 A; 2 (A; 2 + A: 2 ) = IChsin 2 9 cos 9^ in agreement 
with (78) for n = 1. The inversion antisymmetry implies that Z is purely 
imaginary. 

The axisymmetric TSH belonging to the irreducible representation of di- 
mension 9 is 

— 4A; 2 + 3A; 2 + k y %k x k y —8k x k z 

A 9 = 2k x k y — \k\ + 3ky + A; 2 —8k y k z (83) 

—8k x k z -8k y k z 8k 2 z - 4 (fc 2 + A; 2 ) _ 

so that 

A 9 : P = k~ A (24A: 2 (A; 2 + A; 2 ) - 3(A; 2 + k 2 y ) 2 - 8kf) (84) 

and 

Z = AT 4 (30A; 2 (A; 2 + k 2 ) — 5(A; 2 + k 2 ) 2 ^j = 5 sin 2 9(7 cos 2 9 — 1) (85) 

in agreement with (78) provided terms with both n = 0 and n = 2 are 
included. 

8 Conclusions and perspectives 

We summarize the main points of this paper as follows: 
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1. The description of weak anisotropy by tensor functions of wavenumber 
k alone requires two tensors as in Eq. (38). These tensors represent the 
distinct effects of directional anisotropy and polarization anisotropy as 
defined by Cambon et al. [8, 9, 10]. The reduction to a single ten- 
sor function of k compromises the kinematics by introducing implicit 
assumptions about anisotropy, usually that directional anisotropy van- 
ishes. We emphasize in particular that the correlation tensor cannot be 
uniquely and self-consistently reconstructed in terms of the Reynolds 
stresses alone. 

2. Some special cases of anisotropy were considered: short time response 
of turbulence to arbitrary strain and the long time nonlinear response 
of turbulence to small strain. Polarization and directional anisotropy 
are related differently in each limit; consequently, a model based on a 
single spherical average cannot be consistent with both limits. Turbu- 
lence under rapid rotation leads to a conclusion opposite to that usually 
adopted in models, namely dominant directional anisotropy and van- 
ishing polarization. 

3. The description of anisotropy by two tensor functions Y\^ e,z \k,t) was 
exhibited as the lowest order in an infinite expansion in scalar and 
tensor spherical harmonics generated by the SO 3 decomposition. In 
this case, polarization and directional anisotropy are both described 
by a traceless second rank tensor, but the description becomes more 
complex at higher order. 

To conclude, we note some important open issues for anisotropic turbu- 
lence. The first is the question of connecting the higher order coefficients 
in SSH and TSH expansions to spherically averaged higher order moments 
of the spectral tensor. A related question is whether it is possible to model 
the higher order terms by tensor products of H^, H^ z \ and the mean ve- 
locity gradient dUi/dxj if a mean flow is present. Recent studies in rotating 
(Bellet et al. [21]) and/or in stably stratified turbulence have confirmed that 
the anisotropy identified by the angle distribution of £{k ) can be very large, 
and is only reflected by small or moderate values of H-j or In such 
cases, the SSH expansion of R would need to be carried out to extremely 
high order. This expansion is unlikely to be practical; this raises the ques- 
tion whether a different characterization of large anisotropy than spherical 
harmonics expansions may be required. 
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